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ABSTRACT 

We generate random Gaussian turbulent velocity fields with a Kolmogorov spec- 
trum and use these to obtain synthetic line-of-sight velocity profiles. The profiles are 
found to be similar to line profiles observed in molecular clouds. We suggest methods 
l/^ for analysing measured line profiles to test whether they might arise from Gaussian 

Kolmogorov turbulence. 
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1 Introduction 

O Maps of the velocity fields of molecular clouds (Falgarone & Phillips 1990) present a pic- 

^ ture which is reminiscent of turbulence seen in fluid simulations (Porter, Pouqet & Woodward 

l/^ 1994; Falgarone et al. 1994). Velocity variations Av are found to depend on scale A following 

a self-similar relation, Av oc A", where a lies between 0.3 and 0.5. A particularly simple 
picture of turbulence in a homogeneous, incompressible fluid was described by Kolmogorov 
I (1941). In this picture, a fluid which is agitated by some process on a given scale sets up a 

cascade of turbulent vortices of diminishing size, where the characteristic velocity on scale A 
c/3 goes as vx ~ A^/"^. In the language of spectral analysis, the Kolmogorov velocity field corre- 

^ sponds to a Gaussian random field with a steep power spectrum, Pv{k) = ~ k~^^^^. 

Modern theories of turbulence have gone much beyond Kolmogorov's ideas, with increasing 
emphasis being placed these days on phenomena related to "intermittency." Nevertheless, 
we feel that it is useful to compare the turbulence characteristics of molecular clouds to Kol- 
mogorov's original ideas since this is still the simplest description that we have of a turbulent 
fiuid. 

In this paper, we are interested in testing the possibility that the observed velocity 
fields of molecular clouds follow a Kolmogorov spectrum. To do this we set up 3 dimensional 
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random realizations of turbulent velocity fields and then synthesize observed maps of velocity 
profiles by projecting the data down one face of the cube assuming that the fiuid is of 
constant density and optically thin. We use these maps of line profiles to demonstrate new 
observational tests for the Kolmogorov hypothesis. We also compare our results to the similar 
analysis by Falgarone et al (1994) of the 3D fiuid simulations of compressible turbulence of 
Porter et al. (1994). We defer a full analysis of molecular cloud data to a future paper. 

2 A Realization of a Turbulent Field 
2.1 Setting up the Field 

Our goal is to set up a random realization of a turbulent velocity field. The Kolmogorov 
power spectrum of the velocity field for homogeneous incompressible turbulence is 

p^k) = {\vk\') ~ (1) 

where k is the wave number. Since the fiuid is incompressible, we recognize that the fiow is 
divergence free, V • v = 0. Therefore we may define a vector potential. A, from which we 
can derive the velocity field through v = V X A. The components of A are again described 
by a Gaussian random field. By dimensional arguments, the corresponding power spectrum 
is 

ilM') k-''/\ (2) 

The dispersion of |A| at a field point diverges for a power law this steep, so we impose a 
small k cut-off to assure convergence. We introduce a cut-off wavenumber k^in and redefine 
the vector potential power spectrum as 

{\A,\') = C{e + kiJ-'y\ (3) 

where C is a constant. The scale Xmax = '^'^/kmin represents the "outer scale" of the 
turbulence, the scale on which the medium is stirred by whatever feeds the turbulence. 
In an interstellar molecular cloud we expect Xmax to be comparable to the size of the cloud. 
In our simulations, we take Xmax equal to the size of the computational box. 

We generate random realizations of the turbulent velocity field by borrowing the meth- 
ods used in numerical cosmology to set up Gaussian random density fields (e.g. Efstathiou et 
al. 1985). We randomly sample components of the Fourier transform of the vector potential, 
Afc, in a periodic cube with equally spaced lattice points. The real and imaginary compo- 
nents of Afc at coordinate [k^^ ky, k^) are sampled by selecting an amplitude from a Rayleigh 
distribution with dispersion given by the power spectrum in equation (3) and a phase angle 
uniformly distributed between and 27r. (This is equivalent to sampling from a cylindrical 
bivariate Gaussian.) The Fourier transform of the velocity field is then computed by 
taking the curl of A in the Fourier domain: 

Vfc = -ik X Afc. (4) 

In practice, we only need the line-of-sight component of the velocity field for comparison to 
observations so only two components of the vector potential are generated. The velocity field 
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in real space is finally obtained by taking the inverse transform of v^. We do the transform 
using a fast Eourier transform (EET) inside a 256"^ box. 

2.2 Simulated Observations 

Once we have generated the velocity field, we make synthetic velocity profiles by binning 
the velocities with equal weighting in X X 256 unit volumes across the face of the cube, 
as described by Ealgarone et al. (1994). We have used = 8 and 16 to produce 32 X 32 
and 16 X 16 pixel maps (Eigure 1). These velocity profiles would correspond to measured 
line profiles for a homogeneous cloud that is optically thin. 

Since our velocity field has been generated through a Gaussian random process, the 
average of all the velocity profiles should be a Gaussian with dispersion, 

""' = 2^/0 ^'^^^-(^) = i0vri/2r(2/3)r(5/6)^"™'- 

Eor the realizations, we have used C = 1 and kmin = 2Tr / L where L is the boxlength set to 
L = 1. The theoretical and measured dispersions of the line-of-sight velocity agree to 1% and 
are (Ty = 0.13. The velocity profiles in Eigure 1 are just probability distribution functions 
(PDEs) of the velocity along the line-of-sight. The map shows that many of the velocity 
profiles have strongly non-Gaussian features such as skewness, kurtosis, double peaks etc. 
These velocity maps closely resemble the intermediate time results of the fiuid simulations of 
Ealgarone et al., confirming that they were seeing the development of a Kolmogorov turbulent 
cascade in their simulations. Our results also visually resemble the maps of real molecular 
clouds shown in Ealgarone & Phillips (1990). Because of the steepness of the Kolmogorov 
spectrum, our purely random realizations are able to produce quite a variety of distortions of 
the line profiles just as in the observations. Based on this, it would appear that the observed 
profiles are qualitatively similar to a Gaussian random turbulent velocity field. 

Not only do the velocity profiles in Eigure 1 deviate strongly from a Gaussian but there 
is also a spatial correlation of the profile distortions. We have quantified this effect by fitting 
the profiles with an expansion in Gaussian-Hermite (GH) polynomials. These functions are 
an orthonormal basis set which are usually applied in the context of the harmonic oscillator 
in quantum mechanics, but they are also handy fitting functions for our purpose. There are 
different definitions for the GH polynomials but we use the definition prescribed by Van der 
Marel & Eranx (1993) in their analysis of the line profiles in galaxy spectra. The line profile 
is first fit to a Gaussian to determine the best fit of the mean velocity, U, and the dispersion, 
cr. With these parameters fixed, the profile is then fit with a higher order expansion, 

/(„, = 2!!WL^f.,^^^„„a (6) 
I ,=3 J " 

where 7 is a normalizing constant, a[w) = (27r(T^)~^/^ exp(— 1/;^/2) is the best-fitting Gaus- 
sian, and Hj[w) are the GH polynomials. The additional parameters we wish to fit are 
the coefficients hj. The coefficients hi and h2 are set to zero because the process of first 
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fig-mapl6.ps 



Fig. 1. -Simulated 16 X 16 map of velocity profiles created by binning velocities in 16 X 16 X 256 
rods on the face of a 256"^ cube. Although the mean of all these profiles is Gaussian, the individual 
profiles show noticeable deviations from the mean Gaussian. 
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Table 1 
Map Coefficients - a = 0.33 



Map 


V 


cr 




h3 




/i4 






mean rms 


mean 


rms 


mean 


rms 


mean 


rms 


16 X 16 


0.000 0.062 


0.111 


0.021 


0.001 


0.054 


-0.008 


0.033 


32 X 32 


0.000 0.070 


0.121 


0.024 


0.001 


0.077 


0.011 


0.053 



fitting tlie Gaussian a[w) to tlie individual profiles includes their effect on the function. In 
general, the even coefficients broaden or narrow the profiles while the odd coefficients add 
skew in either direction depending on the sign. The coefficients and h4 are similar to the 
third and fourth moments of the distribution which Ealgarone et al. used to characterize the 
data. Eigure 2 shows a subsample of the profiles in Eig. 1 fitted up to the fourth order (i.e. 
including and h4). We see that the fits are fairly good. The GH expansion is therefore a 
useful way to characterize the shapes of the line profiles. 

Eigure 3 shows maps of U, cr, and h4 for one of our random velocity fields. The 
velocities are binned in rods of dimension 8 X 8 X 256 elements across the face of the cube to 
produce 32 X 32 velocity profile maps. Maps at 16 X 16 resolution show the same features. 
The map of V immediately shows the effect of imposing a small wave number cut-off, kmin- 
Eor this realization, we selected kmin equal to the maximum wavelength of the box i.e. 
kfnin = 27r/iv, whcrc L is the box length. The steepness of the spectrum suggests that the 
velocity field should be described by a dominant wave of this length with superimposed small 
scale ripples. The V map does show a dominant wave. The portions of the map with mean 
velocities moving away and towards us are spatially connected revealing the dominance of 
one large scale as imposed by kmin- Indeed the net impression is that the fiuid is rotating 
around an axis pointed towards the northeast. This makes one wonder whether some of the 
cloud rotations that have been claimed are merely random turbulent eddies on the outer 
scale of the cloud. The map of a shows that the regions with broader or narrower than 
average peaks tend to be spatially clustered in different regions of the map. The maps of 
and h4 are much noisier and do not show much spatial correlation. These higher order 
effects can apparently vary significantly from point to point. 

Another useful statistic is the map average and rms values of U, a and the coefficients 
and h4 for an ensemble of realizations. This should give a quantitative feel for the degree of 
non-Gaussian deviations to be expected in a map and may be a useful statistic for comparison 
to real data and fiuid simulations. The average of the mean and rms values of the coefficients 
for an ensemble of 10 maps are given in Table 1. The mean velocity is zero as expected and 
the rms of the mean velocity is about half the 3D value of the expected dispersion, (Ty = 0.13. 
The mean value of a for the projected profiles is approximately equal to the 3D value. The 
and h4 map values suggest that there is an approximately 5% to 8% deviation from the 
best-fitting Gaussian. The 16 X 16 map has slightly smaller values for all of the coefficients 
showing the effect of a large effective smoothing radius. 
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fig-map4.ps 



Fig. 2. -A blow-up of the profiles (solid ) in the lower right corner of Fig. 1 with GH expansion fits 
to 4th order (dashed ). The fits are generally good (except for the double-humped profile at the 
bottom), demonstrating that this technique is a good way to quantify the non- Gaussian deviations 
of the line profiles. Higher order expansions can improve the fits but the the statistics of these fits 
should be interpreted with caution. 
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fig-contmap.ps 



Fig. 3. -Contour maps of V , a, h^, and /i4 for a Kolmogorov turbulent velocity field. The contours 
in V extend between -0.20 and 0.20 in steps of .02, a between and 0.30 in steps of 0.02, 
between -0.5 and 0.5 in steps of 0.05 and /i4 between -0.25 and 0.25 in steps of .025. The V map 
shows a dominant wave of length Xmax = '^'^/kmim as expected for this steep power spectrum, with 
regions approaching and receding being spatially connected. The a map shows peaks and valleys 
corresponding to the broadest and narrowest distributions in Fig. 1. The and /i4 maps show 
much less spatial correlation on large scales than V or a. 
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A comparison of the maps of V and can also be used to discriminate between clouds 
that are turbulent and those that are rotating. For a rotating cloud, velocity profiles are 
asymmetric because of an extended low- velocity tail from gas observed in front of and behind 
the tangent point along the line-of-sight. One therefore expects V and to have opposite 
signs from point to point in a rotating cloud; the maps of V and should be similar and 
should resemble each other except for a sign change. Figure 3 shows that when there is 
turbulence the field of is incoherent and does not resemble the more regular V field. 
This diagnostic can in principle be used to distinguish a turbulent medium from one with a 
dominant net rotation. (We thank the referee, R. van der Marel, for pointing this out.) 

Fits with higher order coefficients, hi^ are also possible and may contain additional 
information. We have made fits with coefficients up to hg and have succeeded in fitting the 
double humped profiles in Figure 2. The fitted values of and h4 do not change significantly 
when including these higher order terms so the results of Table I remain unchanged. However, 
care should be taken in interpreting the maps with and higher. A simple mean and rms 
of the map can be misleading, since the distribution of higher order hi across the map can 
be highly non-Gaussian. In most cases where the profiles are featureless and symmetric, the 
values of higher order hi quickly die off. Profiles with many bumps and wiggles tend to have 
large amplitude hi which alternate in sign between even and odd coeffiecients as one goes 
to higher order. The higher order coefficients are therefore not very informative. Fits to h4 
and perhaps as high as hg are all that is required to quantify the non- Gaussian behavior of 
the velocity profiles. 



2.3 A Steeper Spectrum 

Observations of real molecular clouds suggest velocity scaling relations v oc X" with a as 
large as 0.5 (Falgarone & Phillips 1990). For a = 0.5, the velocity power spectrum is steeper, 
going as Pv{k) ~ k~'^. This steeper spectrum should show more non-Gaussian features in the 
line-of-sight velocity profiles than we have seen with the Kolmogorov spectrum. We repeated 
the previous calculations using a k~'^ velocity spectrum, again introducing a vector potential 
with cut-off at kmin as in equation 3. In this case, the rms variation in the velocity is (cf. 
equation 5) 

OTT 

With C=I and k^in = Stt, (Ty = 0.079. Table 2 shows the values of mean and rms map 
coefficients for the steeper spectrum for comparison to Table I. The rms value of V is larger 
relative to a than in the previous case. The rms values of hs and h4 are also larger, indicating 
more non-Gaussian features in the line profiles than expected for the Kolmogorov spectrum. 



3 Comparison with Simulations of Compressible Turbulence 



Falgarone et al. (1994) have recently analyzed Porter et al.'s (1994) three-dimensional 
fiuid simulations of compressible turbulence with a moderate Mach number (initial rms 
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Table 2 





Map Coefficients 


— a = 


0.5 






Map 


V a 


h 


3 








mean rms mean rms 


mean 


rms 


mean 


rms 


16 X 16 


0.000 0.043 0.069 0.013 


0.000 


0.084 


-0.016 


0.059 


32 X 32 


0.000 0.045 0.068 0.015 


0.001 


0.107 


-0.017 


0.081 



value of 1.1) to compare the resulting simulated velocity field to the observations of velocity 
profiles in molecular clouds. They inject energy into the fiuid modelled as a perfect gas 
with 7 = 4/3 at t = and let it evolve freely, following the development of shocks using 
the Piecewise-Parabolic method. A turbulent velocity field quickly develops and energy is 
dissipated in the shocks. At various times, they produce simulated line-of-sight observations 
of their computational box for comparision to molecular cloud observations. 

The velocity profiles are initially non-Gaussian as the gas begins to dissipate after its 
initial agitation. As the gas settles down, the non-Gaussian features gradually disappear. 
The intermediate snapshot in their simulation produces velocity profiles very similar to the 
observations and the power spectrum is dominated by energy in incompressible modes closely 
following a Kolmogorov scaling, E{}i) oc k~^/^ . The close resemblance of their measured 
velocity profiles to those in Eigure 1 further supports the idea that their simulations of 
turbulence are settling to the Kolmogorov spectrum. 

Ealgarone et al. (1994) also pointed out that there was evidence for the process of 
intermittency in their simulations. They measured a non-Gaussian tail in the distribution 
of vorticity, u = |V X v|, as measured along the line of sight as well as in phase space 
plots of the vorticity versus the 3D bulk velocity at different points along the line of sight. 
In particular, the appearance of dense striations extending to large vorticity for fixed bulk 
velocity suggested the localized process of intermittency. They noted that this behavior could 
also be responsible for some of the non- Gaussian features of the velocity profiles though the 
large variations in bulk fiow along the line of sight due to steepness of the Kolmogorov 
spectrum were the dominant effect in producing the profiles. 

We do not expect to see streaks in the u — v phase space with a Kolmogorov spectrum 
since the underlying random process generating the field is Gaussian (though "non-Gaussian" 
velocity profiles are created because of the steep power spectrum). Eigure 4 shows a plot of 
the vorticity, cu, versus the three components of velocity along 3 adjacent lines of sight down 
the X face of the cube. They are chosen because of their pronounced non-Gaussian velocity 
profiles. Also, these synthesized profiles were generated at 128"^ resolution because the extra 
memory overhead to generate, cu, at 256"^ resolution was too great for our local resources. 
These plots are similar to Eigure 9 in Ealgarone et al. (1994). The Kolmogorov turbulent 
fields show no signs of the non-Gaussian streaks in u — v phase space that are seen in the 
fiuid simulations. 

Eigure 5 emphasizes again why we expect non- Gaussian profiles from an underlying 
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fig-rotv.ps 



Fig. 4. -Scatter plot of |rot v\ vs. the velocity components v^, Vy, and along 3 adjacent lines 
of sight averaged including 8x8 pixels on the face of 128"^ random velocity field. Also shown are 
the line-of-sight velocity profiles of the component. There are no signs of the striations in phase 
space seen by Falgarone et al. 1994 in their fiuid simulations which they suggested were signs of the 
non-Gaussian effects of intermittency. The purely Gaussian turbulent field has a bland appearance 
in this representation of phase space. 



TURBULENCE 



11 



Gaussian process. This figure shows the distribution of the iine-of-sight component of veloc- 
ities in the 8x8 pixel region for the three profiles in Eigure 4. The large scale wave of the 
bulk fiow imprinted by the dominant long waves of the Kolmogorov spectrum is once again 
demonstrated. 

4 Conclusions 

Despite the simplifying assumptions of homogeneity, isotropy and incompressibility, 
Kolmogorov turbulence seems to be a good description of the velocity profiles of molecular 
clouds. Erom our analysis, we draw the following conclusions: 

1. The line-of-sight velocity profiles derived from an optically thin Kolmogorov turbulent 
velocity field exhibit non- Gaussian features which resemble the features seen in real 
molecular clouds and simulations of compressible turbulence. These features occur 
because of the steepness of the Kolmogorov power spectrum i.e. a few large scale 
waves dominate the character of the bulk fiow. 

2. The deviation from a Gaussian of the velocity profiles as quantified by higher order fits 
with Gaussian-Hermite polynomials are at the 5% to 10% level (rms value of the 
and /i4 coefficients over the simulated map). A steeper velocity spectrum [Py ~ k~'^) 
creates more non- Gaussian features in the line profiles which show up as larger values 
of and 

3. The presence of an outer scale, Xmax^ to the turbulent field reveals itself as a large 
scale wave in the mean projected velocity, mimicing a net rotation. A comparison of 
the V and fields derived from the line-of-sight profiles can be used to discriminate 
between rotation and turbulence in the cloud: a strong anti-correlation between V and 

indicates that rotation dominates, while a regular V field with an incoherent field 
indicates that the apparent rotation is merely part of a turbulent cascade. 

4. Eluid simulations show evidence of intermittency as revealed by non-Gaussian stripes 
in plots of vorticity versus bulk velocity. These features do not occur in our simple 
model of Kolmogorov turbulence because it is a fundamentally Gaussian model. Given 
that our model produces line profiles that are fairly similar to actual line profiles of 
molecular clouds, it is not clear if the data provide evidence for intermittency. 



We thank the referee R. van der Marel for useful comments and suggestions. TP and 
RN acknowledge NSE grants AST 93-13929 and AST 9148279 respectively. JD recognizes a 
CfA postdoctoral fellowship. 
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fig-vline.ps 



Fig. 5. -The distribution of the line-of-sight velocities down the line- of- sight. The same velocity 
profiles from Figure 4 are shown. The outer scale, Xmax, is apparent in the variation of the bulk 
velocity down the line of sight and is responsible for the observed non-Gaussian velocity profile. 
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